Fourier Analysis of the Parametric Resonance of the Neutrino 
Oscillation in the Presence of Inhomogeneous Matter 



oo 

o 
o 

(N 

> 

o 
O 



in: 



> 
o 

m 
d 

oo 
o 

> 

• 1—1 

x 



Joe Sato a * Masafumi Koike a j Toshihiko Ota b J and Masako Saito a § 

™ Physics Department, Saitama University, 255 Shimo-Okubo, Sakura-ku, Saitama 338-8570, Japan 
b Institiit fur Theoretische Physik und Astrophysik, Universitat Wurzburg, 97074 Wurzburg, Germany 



Abstract 

We study the parametric resonance of the neutrino oscillation through the matter whose density varies 
spatially. The Fourier analysis of the matter effect enables us to clarify the parametric resonance condition, 
which is summarized in a frequency matching between the neutrino oscillation and the spatial variation 
of the matter density. As a result, the n-th Fourier mode of a matter density profile modifies the energy 
spectrum of the — > v c appearance probability at around the n-th dip. 



The oscillation probability of neutrinos passing through the interior of the Earth is strongly affected by the 
matter on the baseline. In such cases, the spatial variation of the matter density can make the probability quite 
different from that under the constant matter density [1, 2, 3]. The possibility of the parametric resonance of 
the neutrino oscillation with the matter density has been discussed [1]. We exploit the Fourier analysis of the 
matter density profile to study the effect of parametric resonance [3, 4]. Generally, the parametric resonance 
of an oscillation emerges when an oscillation parameter changes with a resonant frequency, which is typically 
twice the natural frequency of the system. Our formulation based on the Fourier analysis gives a clear view of 
this frequency matching, which is essential to the parametric resonance. 

We carry out our analysis in a two-flavor framework. This simplification does not spoil our point as can be 
seen in Fig. 1, in which we compare the — > v e appearance probability in the three- and two-flavor calculation, 
taking the baseline length L = 12000 km. The three-flavor result assumes the density profile of the Preliminary 
Reference Earth Model (solid curve), while two-flavor result uses the constant density case (dotted curve) and 
the cosine variation around the average density (broken curve). Renormalized by the factor of sin 2 26*23, the 
two-flavor result with the cosine profile shows the similar feature as the three-flavor one. 

The evolution equation of two flavors of neutrinos v c and in the matter density p{x) is 
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where dm 2 , 0, and E are the mass-square difference, the mixing angle, and the energy of neutrinos, respectively. 
The matter effect a(x) is given by a(x) = 2v / 2 E 'G 'f N xY c p{x) , where G-p is the Fermi constant, Na is the 
Avogadro number, and Y c is the proton-to-nuclcon ratio. Defining A = 6m 2 L/2E, A m (x) = a(x)L/2E, and 
z{x) = u e (x)exp[J dsiA m (s)/2], and introducing an Fourier expansion A m (£) = J^^g A m „ cos 2n7r^ with 
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Figure 1: Oscillation probability from to v c for the baseline length of 12000 km. The solid curve (A) is 
evaluated by three-flavor analysis with the matter density profile taken from the PREM. The dotted (B) and 
broken (C) curves are results of the two flavor calculation with (B) constant density with po = 7.58 g/cm 3 , 
and (C) cosine profile p(x) = pa + pi cos2ttx/L with pi = —4.32 g/cm 3 . Values of other parameters taken are 
described in the figure. 
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Figure 2: The — > z/ e oscillation probabilities with p n = 0, —0.5, and —1.0 g/cm 2 , where (A) n = 1 and (B) 
n = 2. The former includes pi = —5.0 g/cm 2 case also. Values of other parameters are shown in the graphs. 



£ = x/L, we obtain from Eq. (1) 

z" (£) + - [(A m0 - A cos 20) 2 + A 2 sin 2 20 + 2(A m0 - A cos 26) ^ A m „ cos 2n< 
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( y A m „ cos 2n7r£ ) + 4n7ri > A m „ sin 2n7Tc; z(£) = . 
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The parametric resonance can arise under the presence of the density profile. Suppose that only the n-th 
Fourier mode of the density profile is present so that A m (£) = A m o + A m „ cos2n7Tc;. The parametric resonance 
is expected when the frequency of the density profile 2nir is twice the natural frequency of the system luq = 



(A m o — A cos 2d) 2 + A 2 sin 2 26 — h YlnLi A 2 m /2 , i.e. loq = 2rnr/2 = mr. This frequency matching condition 
amounts to 

E = E^= 6 -^ . (3) 

2 A m0 cos 26 ± yV^ 2 - A 2 n0 sin 2 26 - \ EZi 

We show in Fig. 2 the — > v c oscillation probability for various values (A) of pi and (B) of p^. Figure 2 
indicates that the n-th Fourier mode of the matter density modifies the oscillation probability around the n-th 
dip: u>o = n7r, which corresponds to the parametric resonance energy Eq. (3). This modification is understood 
by the parametric resonance as shown elsewhere [5]. 
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In summary, wc studied the parametric resonance of the neutrino oscillation through the matter whose 
density varies spatially, analyzing the matter effect by the Fourier expansion. We have shown that the n-th 
Fourier mode of the density profile modifies the oscillation probability around the parametric resonance energy 
defined in Eq. (3). This condition is understood in terms of the frequency matching between the natural 
frequency of the system and that of the density profile. 
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